Numerical evidence for 'multi-scalar stars' 
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We present a class of general relativistic soliton-like solutions composed of multiple minimally 
coupled, massive, real scalar fields which interact only through the gravitational field. We describe 
a two-parameter family of solutions we call "phase-shifted boson stars" (parameterized by central 
density po and phase 5) , which are obtained by solving the ordinary differential equations associated 
with boson stars and then altering the phase between the real and imaginary parts of the field. These 
solutions are similar to boson stars as well as the oscillating soliton stars found by Seidel and Suen 
[E. Seidel and W. M. Suen, Phys. Rev. Lett. 66, 1659 (1991)]; in particular, long-time numerical 
evolutions suggest that phase-shifted boson stars are stable. Our results indicate that scalar soliton- 
like solutions are perhaps more generic than has been previously thought. 
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I. INTRODUCTION 

The nature of dark matter in the universe is currently 
an open question in physics, with many models being 
proposed to fill this gap in our understanding, some of 
which resort to the use of exotic matter. Of interest to us 
is one class of models composed of massive scalar fields 
coupled to the general relativistic gravitational field, from 
which compact, star-like solutions can be formed, solu- 
tions which go by the names of "oscillating soliton stars" 
(or "oscillatons" ) [H |^ for real fields and "boson stars" 
||, |, |, I, 0, D for complex fields. These star-like so- 
lutions have received renewed attention recently, and a 
substantial body of evidence has been advanced in an ef- 
fort to show that these fields may be key players on both 
galactic |l^ and cosmological scales. Boson 

stars have been suggested as alternatives to primordial 
black holes as well as supermassive black holes in 
the centers of galaxies , and their gravitational lens- 
ing properties have been explored jl5| , further developing 
the treatment of these solutions as objects of astrophys- 
ical interest. 

Apart from the possible astrophysical relevance of 
these star-like objects, we find them interesting to study 
from a mathematical standpoint as well, for their proper- 
ties as soliton-like solutions in the nonlinear dynamics of 
general relativity. The "solution space" of general rela- 
tivity is still largely unexplored, and these scalar objects 
comprise simple systems with which to conduct investi- 
gations. It is from this viewpoint that we will proceed in 
this paper. 

In 1991, Seidel and Suen ^ considered the model of a 
real massive scalar field, minimally coupled to the gen- 
eral relativistic gravitational field, with the additional 
simplifying assumption of spherical symmetry. These au- 



thors were interested in the existence of "nontopological 
solitons" in the model: that is, whether the equations 
of motion admitted stable, localized, non-singular distri- 
butions of matter which could be interpreted as "scalar 
stars". A theorem due to Rosen suggested that, 
should such solutions exist, they could not be static. 
Thus, Seidel and Suen looked for periodic configurations 
by substituting a particular Fourier ansatz into the equa- 
tions of motion, and solving the resulting hierarchy of 
ODEs via a generalized shooting technique. The authors 
found strong evidence that periodic star-like solutions did 
exist, and, via direct numerical simulation, demonstrated 
that their "oscillating soliton stars" , if not absolutely sta- 
ble, had lifetimes many orders of magnitude longer than 
the stars' intrinsic dynamical times. ^ 

These results were surprising to some researchers, par- 
ticularly since the model has no conserved Noether cur- 
rent, which, it had been argued, was responsible for the 
existence of solitonic solutions in other non-linear field 
theories involving scalar fields ||l^, |l8). However, at least 
heuristically, we can understand the existence of the os- 
cillating stars as arising from a balance between the at- 
tractive gravitational interaction and the effective repul- 
sive self-interaction generated by the mass of the scalar 
field {i.e. via the dispersion relation of the Klein-Gordon 
equation) . 

Recently it was shown by Ureha-Lopez ||] that ap- 
proximate solutions for boson stars and oscillating soli- 
ton stars, or "oscillatons" as he calls them, can both be 
derived from a single set of equations in a sort of "sta- 
tionary limit" . The similarities seen between boson stars 



*E-mail address: shawley@aei.mpg.de 
tE-mail address: choptuik @ physics. ubc.ca 



^ More precisely, the oscillating stars constitute a one-parameter 
family which may be parametrized by the mean (period- 
averaged) central density, po- As with other relativistic stellar 
models, a plot of total (ADM) mass versus po exhibits a maxi- 
mum at pg , which seems to coincide with a transition from stable 
to unstable configurations. As expected, only stars with po < Pg 
could be stably evolved for long times. 
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and oscillating soliton stars in terms of their curves relat- 
ing mass, radius and central density can thus be related 
formally. 

In this paper, we build on the works of Seidel and Suen 
and Urefia-Lopez by considering a matter content con- 
sisting of multiple scalar fields, and we find some further 
ways in which boson stars and oscillating soliton stars 
are similar. For the specific case of two scalar fields, we 
find evidence for a new family of quasi-periodic, solitonic 
configurations.^ Together with previous results, this sug- 
gests that solitonic solutions are generic to models which 
couple massive scalar fields through the Einstein gravita- 
tional field. We should note, however, that we make no 
attempt to address the question of whether these mul- 
tiple scalar fields actually exist in nature; rather we are 
interested in their existence as valid mathematical solu- 
tions in the Einstein-Klein-Gordon system. 

We begin by considering n real, massive Klein- 
Gordon fields i = 1,2, .. .n, without additional self- 
interaction, minimally coupled to the general relativistic 
gravitational field. Specifically, choosing units such that 
c = 1 and G = 1 , the Lagrangian density for the coupled 
system is 



C = 



R 



E 

i=l 



(1) 



where g = det^^jy, R is the Ricci scalar, and nii is the 
mass of the z-th scalar field. 

We now restrict attention to spherical symmetry, and 
adopt the "polar/areal" coordinate system, so that the 
metric takes the form: 



-a^{t, r) + a^{t, r) dr^ + dil^ 



(2) 



The complete evolution of the scalar fields and spacetime 
can then be given in terms of a Klein-Gordon equation for 
each of the , and two constraints derived from the Ein- 
stein field equations and the coordinate conditions used 
to maintain the metric in the form (||). We solve these 
equations using the same scheme adopted for the criti- 
cal phenomena study described in ]l9t , and only briefly 
review that scheme here. 

We define the following auxiliary scalar field variables: 



a 



(3) 



^ Although the oscillations appear to be periodic, we do not have a 
proof of their periodicity and so we use the term "quasi-periodic" 
to describe their temporal behavior. Small departures from strict 
periodicity over long times scales are visible in the simulation 
results; however, these departures become smaller as we decrease 
the radial mesh spacing Ar, and thus it is plausible that in the 
limit Ar — > 0, the solutions are truly periodic. 



where all variables are functions of t and r, ' = d/dt and 
' = d/dr. The Klein-Gordon equation is written as the 
following system: 

lit = — [ III I - mfaa(j)i 



/ (4) 



where r = rmax is the outer boundary of the computa- 
tional domain. The constraint equations are the "Hamil- 
tonian constraint" 



2r 



i=l 



(5) 



and the "slicing condition" 



a2 - 1 a' 



a — a 



(6) 



For diagnostic purposes, we have also found it useful to 
compute and monitor the quantities, Mi(t, '"max), defined 
by 



Jo 



where 
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2o? 



(8) 



Loosely speaking, we can interpret r,„ax) as the 

total contribution of field i to the ADM mass of the 
spacetime. In particular, so long as no matter out-fluxes 
through r = r,„ax: we have 



const. 



(9) 



We solve Eqs. (^)-(|^) subject to the the boundary 
conditions a{t, 0) = 1 (local flatness at the origin) and 
a{t, Tinax) = ?'max) (so that t mcasurcs proper time 

as r ^ 00). As in |l9| , we use the Sommerfeld condition 
for a massless field to set the values r,„ax), r,„ax) 
and n(i, r,„ax)- Since the Sommerfeld condition is not 
ideal for a massive field, we ran our simulations with dif- 
ferent values of r„iax, testing for any periodicity or other 
effect which might be a function of rmax, and usually ran 



with an 



which was large compared to the time for 



which we ran the simulation. Even with smaller rmax, we 
found our results to be essentially independent of Tmax 
and attribute this to the fact that there is very little 
scalar radiation emitted from the soliton-like objects con- 
sidered here. Our results are also essentially independent 
of the resolution of the finite differencing algorithm and 
the Courant-Friedrichs-Levy factor, At/Ar, and we con- 
firm that our results converge in a second-order-accurate 
manner using independent residual evaluations. 
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II. 



"PHASE-SHIFTED BOSON STARS" 



We start by considering the case n — 1, so that our 
matter content is a single scalar field, (j>{t,r). We note 
that the Hamiltonian constraint (^) and the slicing con- 
dition (j^) are unchanged if we trivially decompose (j) into 
two identical fields {i.e. now choosing n = 2), namely 
(f)i{t,r) and (j)2{t,r) — (j)i{t,r) (with 7712 — mi — 1), such 
that 



1 
71 



+ 02). 



(10) 



Further, we note that for fixed a(i, r) and a(t, r), if 4>{t, r) 
is a solution of the Klein-Gordon equation (Q), then so 
is K(j){t,r) where n is an arbitrary real constant. Since 
a soliton solution of the system (^-(||) is the oscillat- 
ing soliton star, we see that a trivial multi-scalar soliton 
solution can be obtained by constructing an oscillating 
soliton star with a single field (f> (as described in 0]) 
and then reinterpreting it as a two-field solution in which 

Moreover, if we wish to model a boson star |^, |[ ^, 0] 
with no self-interaction potential (often called a "mini- 
boson star" as in jsj), then we have one massive com- 
plex scalar field (j>{t,r), for which the real and imag- 
inary parts behave like two real-valued scalar fields: 
4>{t,r) — (j)i(t,r) + i(j)2{t,r). The boson star ansatz is 
(j){t, r) — 4>{r) exp(±ict;t), where (j){r) is real. This implies 



hit,r) = (^(r) cos([ji) 
h{t,r) = 4>{r) cos{ujt + S) , 



(11) 



where 6 = '^-k 12. 

Thus we see that soliton stars and boson stars can both 
be obtained using two real scalar fields with a constant 
temporal phase difference. For soliton stars, the fields 
are identical for all r and i; whereas for boson stars, the 
fields have identical r-dependence, and the ^-dependence 
is the same to within a phase. (In each case, the central 
density po is uniquely fixed by the value of the field at 
the origin, e.g. 0(0).) 

The work described in the remainder of this paper be- 
gan in the midst of our numerical evolutions of boson 
stars [|l9|. The question arose, "What happens if we solve 
for the boson star initial data, then 'manually' change 
the phase between the two fields, and then re-solve the 
constraints to obtain the metric variables?" For future 
reference, we term such a modified-boson-star configu- 
ration a "phase-shifted boson star." This modification 
to the boson star data was motivated more by "practi- 
cal" reasons that "physical" ones — we were interested in 
studying oscillating soliton stars but found them difficult 
to construct. 

Taking the boson star initial data and setting 
(/>2(0,r) = (/)i(0, r) resulted in what might be termed a 
"poor man's soliton star" (PMSS). In Figures |l|, 2 and 3, 
we show that the resulting solution is very similar to the 
true soliton star solution, and can perhaps best be re- 
garded as a soliton star with a small perturbation added. 
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FIG. 1: A comparison between initial data for a true soli- 
ton star (SS, solid line) with the "poor man's soliton star" 
(PMSS, circles), which is a phase-shifted boson star where 
the two fields are in phase (i.e., 5 = 0), for a particular choice 
of the central value of the field 0(0,0). In this figure, we 
compare the scalar field of the soliton star with one of the 
two (identical) fields comprising the PMSS, where we have 
divided the soliton star field by \/2 in keeping with the rela- 
tion (|l|). The relative difference between the two solutions 
is plotted in the inset. Given that these two solutions are 
obtained by solving two rather different sets of ODEs (three 
simple ODEs for the PMSS, and a complicated system of 10 
ODEs for the SS), we find it remarkable that they are so sim- 
ilar. (We note that although we focus on the scalar fields in 
this figure, the metric variables for the PMSS are also close 
to those of the SS.) 



We then took the boson star initial data (^(r) and dis- 
tributed it to 01 and 02 via Eq. (^l|) using some different 
value of (5, such as i5 = — 7r/6. (Note that we only apply 
Eq. ( pT| ) for the initial data, i.e. at f = 0; one cannot ex- 
pect Eq. ( [Til ) to describe the fields for all t if S ^ i'l'/Z, 
as the S — case demonstrates.) One aspect of the evo- 
lution for such a system can be seen in Figure |^. 

For each of the many values of S we tried, we found an 
apparently stable solution which oscillated in some essen- 
tially periodic manner for very long times. (The phase 
was preserved throughout the evolution, i.e. it is not 
the case that the system reverted to a simple "perturbed 
boson star" over time.) These results led us to conjec- 
ture that there may exist a continuous family of periodic 
soliton-like solutions (parameterized by the phase S) of 
which our "phase-shifted boson stars" are perturbations. 
We hope in the future to construct such a family directly 
via a periodic ansatz of the form used by Seidel and Suen 
for their oscillating soliton stars (with additional terms 
incorporated to account for the nonlinear coupling be- 
tween the two scalar fields). 

We would like to mention that, for a given 0(0), vary- 
ing S has very little effect on the metric variable a(0, r). 
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FIG. 2; Central value of the fields <^i(i,0) = 02(t,O) vs. 
time t, for the "poor man's soliton star" (PMSS). In the top 
panel we show both the evolution of the PMSS field (solid 
line) and, for comparison, a similar evolution for a soliton 
star field (dashed line). In the lower panel we show only the 
"envelope" of the oscillations in the PMSS scalar field; the 
period of variations in the envelope is roughly 32 times the 
intrinsic period of field. Stable evolutions of this system for 
t > 20000 have been obtained. 



FIG. 4: Time development of the quantities Mi{t, rmax) (solid 
line) and M2(t,rinax) (dashed line), defined by (Q), for a 
phase-shifted boson star with 5 — —tt/6. We interpret the 
behaviour seen in this figure as the periodic transfer of sig- 
nificant amounts of energy from one scalar field to the other 
via the intermediary of the gravitational field. This is remi- 
niscent of "beats" in weakly-coupled harmonic oscillators; in 
this case, the coupling between the two "oscillators" is gravi- 
tational. 
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FIG. 3: Fourier transform of the PMSS evolution shown in 
Figure ^ The spikes in the spectrum correspond closely to the 
harmonics {lo,3uj,5uj, ...} of the oscillating soliton star. The 
value of u) is an eigenvalue of the soliton star ODE problem; 
for this simulation, u ~ 0.143. 



Consequently the so-called "stability curves" relating to- 
tal mass, radius and central density are essentially the 
same as the curves for boson stars {5 ~ ±tt/2); conse- 
quently, we do not consider it relevant to plot them here. 



III. CONCLUSIONS 

We have strong numerical evidence for the existence 
of a two-parameter family of soliton-like solutions to the 
Einstein-Klein-Gordon system (parameterized by central 
density po and phase 6) for the case of two scalar fields. 
We speculate that extending the system to more scalar 
fields will yield similar results (especially if one uses a 
trivial extension like Eq. ( p^ ). The solutions we refer 
to as "phase-shifted boson stars" consist of boson star 
initial data for which the phase difference S between the 
real and imaginary components of the field have been 
altered. These solutions oscillate in a seemingly periodic 
manner for very long times; thus they appear to be stable. 
For the case of ^ = 0, we obtain close approximations to 
the oscillating soliton stars of Seidel and Suen. For other 
values of 6, we find solutions which also appear to be 
stable and periodic; furthermore we can see mass-energy 
being exchanged between the two fields. 
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